Possible realization of an ideal quantum computer in Josephson junction array 
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We introduce a new class of Josephson arrays which have non-trivial topology and exhibit a novel 
state at low temperatures. This state is characterized by long range order in a two Cooper pair 
condensate and by a discrete topological order parameter. These arrays have degenerate ground 
states with this degeneracy 'protected' from the external perturbations (and noise) by the topological 
order parameter. We show that in ideal conditions the low order effect of the external perturbations 
on this degeneracy is exactly zero and that deviations from ideality lead to only exponentially 
small effects of perturbations. We argue that this system provides a physical implementation of an 
ideal quantum computer with a built in error correction and show that even a small array exhibits 
interesting physical properties such as superconductivity with double charge, 4e, and extremely long 
decoherence times. 



I. INTRODUCTION 



Quantum computing^'^ is in principle a very power- 
ful technique for solving classic 'hard' problems such 
as factorizing large numbers'^ or sorting large lists^. 
The remarkable discovery of quantum error correction 
algorithms^ shows there is no problem of principle in- 
volved in building a functioning quantum computer. 
However, implementation still seems dauntingly difficult: 
the essential ingredient of a quantum computer is a quan- 
tum system with 2^ (with K ^ 104 ~ 106) quantum 
states which are degenerate (or nearly so) in the ab- 
sence of external perturbations and are insensitive to the 
'random' fluctuations which exist in every real system, 
but which may be manipulated by controlled external 
fields with errors less 10^'' — 10^® (big system sizes, K, 
are needed to correct the errors, for smaller errors the 
size of the system, K, gets much smaller)®. Insensitivity 
to random fluctuations means that any coupling to the 
external environment neither induces transitions among 
these 2^ states nor changes the phase of one state with 
respect another. Mathematically, this means that one 
requires a system whose Hilbert space contains a 2^- 
dimensional subspace (called 'the protected subspace'^) 
within which any local operator O has (to a high accu- 
racy) only state-independent diagonal matrix elements: 
{n\0\m) = OoSmn + o(exp(— L)) where L is a parame- 
ter such as the system size that can be made as large 
as desired. It has been very difficult to design a system 
which meets these criteria. Many physical systems (for 
example, spin glasses*) exhibit exponentially many es- 
sentially degenerate states, not connected to each other 
by local operators. However, the requirement that all 
diagonal matrix elements are equal (up to vanishingly 
small terms) is highly nontrivial and puts such systems 
in a completely new class. Parenthetically we note that 
such systems were discussed in philosophical terms by L 
Kant (who termed them noumcnons, in his thinking the 
noumenal world is impenetrable but contains comprehen- 
sible information) in^. 



One very attractive possibility, proposed in an impor- 
tant paper by Kitaev^ and developed further in^" involves 
a protected subspace^^'^^ created by a topological degen- 
eracy of the ground state. Typically such degeneracy 
happens if the system has a conservation law such as the 
conservation of the parity of the number of 'particles' 
along some long contour. Physically, it is clear that two 
states that differ only by the parity of some big number 
that can not be obtained from any local measurement are 
very similar to each other. The model proposed in^ has 
been shown to exhibit many properties of the ideal quan- 
tum computer; however before now no robust and prac- 
tical implementation was known. In a recent paper we 
and others proposed a Josephson junction network which 
is an implementation of a similar model with protected 
degeneracy and which is possible (although difficult) to 
build in the laboratory 

In the present paper we develop ideas of^'^ proposing 
a new Josephson junction network that has a number 
of practical advantages, (i) This network operates in a 
phase regime (i.e. when Josephson energy is larger than 
the charging energy), which reduces undesired effects of 
parasitic stray charges, (ii) All Josephson junctions in 
this array are similar which should simplify the fabrica- 
tion process, (iii) This system has 2^ degenerate ground 
states 'protected' to even higher extent than in^^: matrix 
elements of local operators scale as e^, where e < 0.1 is a 
measure of non-ideality of the system's fabrication (e.g. 
the spread of critical currents of different Josephson junc- 
tions and geomertical areas of different elementary cells 
in the network) (iv) The new array does not require a fine 
tuning of its parameters into a narrow region. The rele- 
vant degrees of freedom of this new array are described 
by the model analogous to the one proposed in^. Even 
when such system is small and contains only a few 'pro- 
tected' states its physical properties are remarkable: it 
is a superconductor with the elementary charge 4e and 
the decoherence time of the protected states can be made 
macroscopic allowing 'echo' experiments. 

Below we first describe the physical array, identify its 
relevant low energy degrees of freedom and the mathe- 
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matical model which describes their dynamics. We then 
show how the protected states appear in this model, de- 
rive the parameters of the model and identify various 
corrections appearing in a real physical system and their 
effects. Finally, we discuss how one can manipulate these 
states in a putative quantum computer and the physical 
properties expected in a small arrays of this type. 

II. ARRAY 

The basic building block of the lattice is a rhombus 
made of four Josephson junctions with each side of the 
rhombus containing one Josephson contact, these rhombi 
form a hexagonal lattice as shown in Fig. 1. We denote 
the centers of the hexagons by letters a,b . . . and the 
individual rhombi by (a6), (cd) . . ., because each rhom- 
bus is one-to-one correspondence with the link (ab) be- 
tween the sites of the triangular lattice dual to the hexag- 
onal lattice. The lattice is placed in a uniform mag- 
netic field so that the flux through each rhombus is $o/2. 
The geometry is chosen in such a way that the flux, 
through each David's star is a half-integer multiple of $0: 
<1> = (ris -I- i)<i>o.'^'' Finally, globally the lattice contains 
a number, K, of big openings (the size of the opening 
is much larger than the lattice constant, a lattice with 
K = 1 is shown in Fig. la). The dimension of the pro- 
tected space will be shown to be equal 2^ . The system 
is characterized by the Josephson energy, Ej = |^/c, of 
each contact and by the capacitance matrix of the islands 
(vertices of the lattice). We shall assume (as is usually 
the case) that the capacitance matrix is dominated by the 
capacitances of individual junctions, we write the charg- 
ing energy as Ec = The 'phase regime' of the net- 
work mentioned above implies that Ej > Ec ■ The whole 
system is described by the Lagrangian 

= TTE^^'i'i ~ + cos{(l)^ - <j>j - fly) (1) 
w) 

where (jji are the phases of individual islands and a.y are 
chosen to produce the correct magnetic fluxes. The La- 
grangian (1) contains only gauge invariant phase differ- 
ences, 4'ij — (pi — (pj — aij, so it will be convenient some- 
times to treat them as independent variables satisfying 
the constraint 4>ij = 27r<I>r/*J'o + 27rn where the sum 
is taken over closed loop F and n is arbitrary integer. 

As will become clear below, it is crucial that the de- 
grees of freedom at the boundary have dynamics identical 
to those in the bulk. To ensure this one needs to add ad- 
ditional superconducting wires and Josephson junctions 
at the boundary. There are a few ways to do this, two 
examples are shown in Fig. la and Fig. lb: type I bound- 
ary (entire length of boundaries in Fig. la, parts AB, 
CD) and type II boundary {BC, AD). For both types of 
boundaries one needs to include in each boundary loop 
the flux which is equal to .Z;, * 7r/2 where Zf, is coordi- 
nation number of the dual triangular lattice site. For in- 
stance, for the four coordinated boundary sites one needs 
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FIG. 1: Examples of the proposed Josephson junction ar- 
ray. Thick lines show superconductive wires, each wire con- 
tains one Josephson junction as shown in detailed view of 
one hexagon. The array is put in magnetic field such that 
the flux through each elementary rhombus and through each 
David's star (inscribed in each hexagon) is half integer. Thin 
lines show the effective bonds formed by the elementary 
rhombi. The Josephson coupling provided by these bonds 
is TT-periodic. a. Array with one opening, generally the effec- 
tive number of qubits, K is equal to the number of openings. 
The choice of boundary condition shown here makes super- 
conducting phase unique along the entire length of the outer 
(inner) boundary, the state of the entire boundary is described 
by a single degree of freedom. The topological order param- 
eter controls the phase difference between inner and outer 
boundaries. Each boundary includes one rhombus to allow 
experiments with flux penetration; magnetic flux through the 
opening is assumed to be + Tn) with any integer m. 

b. With this choice of boundary circuits the phase is unique 
only inside the sectors AB and CD of the boundary; the topo- 
logical degree of freedom controls the difference between the 
phases of these boundaries. This allows a simpler setup of 
the experimental test for the signatures of the ground state 
described in the text, e.g. by a SQUID interference experi- 
ment sketched here that involves a measuring loop with flux 
<l>m and a very weak junction J balancing the array. 



to enclose the integer flux in these contours. In type I 
boundary the entire boundary corresponds to one degree 
of freedom (phase at some point) while type II bound- 
ary includes many rhombi so it contains many degrees of 
freedom. 

Note that each (inner and outer) boundary shown in 
Fig. la contains one rhombus; we included it to allow 
flux to enter and exit through the boundary when it is 
energetically favorable. 



III. GROUND STATE, EXCITATIONS AND 
TOPOLOGICAL ORDER 

In order to identify the relevant degrees of freedom in 
this highly frustrated system we consider first an indi- 
vidual rhombus. As a function of the gauge invariant 
phase difference between the far ends of the rhombus the 
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FIG. 2: The spin degrees of freedom describing the state of the 
elementary rhombi are located on the bonds of the triangular 
lattice (shown in thick lines) while the constraints are defined 
on the sites of this lattice. The dashed line indicates the 
boundary condition imposed by a physical circuitry shown in 
Fig. la. Contours 7 and 7' are used in the construction of 
topological order parameter and excitations. 



potential energy is 

= -2Ej{\ cos(0y/2)| + I sin(</.y/2)|). (2) 

This energy has two equivalent minima, at (jjij = ±7r/2 
which can be used to construct elementary unprotected 
qubit, see-'^^. In each of these states the phase changes 
by ±7r/4 in each junction clockwise around the rhombus. 
We denote these states as | t) and | |) respectively. In 
the limit of large Josephson energy the space of low en- 
ergy states of the full lattice is described by these binary 
degrees of freedom, the set of operators acting on these 
states is given by Pauli matrices cr^}^'^. We now combine 
these rhombi into hexagons forming the lattice shown in 
Fig. 1. This gives another condition: the sum of phase 
differences around the hexagon should equal to the flux, 
$3 through each David's star inscribed in this hexagon. 
The choice — tt/2 is consistent with flux $5 that is 
equal to a half integer number of flux quanta. This state 
minimizes the potential energy (2) of the system. This 
is, however, not the only choice. Although flipping the 
phase of one dimer changes the phase flux around the 
star by tt and thus is prohibited, flipping two, four and 
six rhombi is allowed; generally the low energy configu- 



rations of U{(j)) satisfy the constraint 

b 

where the product runs over all neighbors, b, of site a. 
The number of (classical) states satisfying the constraint 
(3) is still huge: the corresponding configurational en- 
tropy is extensive (proportional to the number of sites). 
We now consider the charging energy of the contacts, 
which results in the quantum dynamics of the system. 
We show that it reduces this degeneracy to a much 
smaller number 2^. The dynamics of the individual 
rhombus is described by a simple Hamiltonian H — ta^ 
but the dynamics of a rhombus embedded in the array 
is different because individual flips are not compatible 
with the constraint (3). The simplest dynamics compat- 
ible with (3) contains flips of three rhombi belonging to 

the elementary triangle, (a,5, c), Qi^abc) — ^Zb'^bc^ca 
therefore the simplest quantum Hamiltonian operating 
on the subspace defined by (3) is 

^ = E Qi-bc) (4) 
[ahc] 

We discuss the derivation of the coefficient r in this 
Hamiltonian and the correction terms and their effects 
below but first we solve the simplified model (3,4) and 
show that its ground state is "protected" in the sense de- 
scribed above and that excitations are separated by the 
gap.'' 

Clearly, it is very important that the constraint is im- 
posed on all sites, including boundaries. Evidently, some 
boundary hexagons are only partially complete but the 
constraint should be still imposed on the corresponding 
sites of the corresponding triangular lattice. This is en- 
sured by additional superconducting wires that close the 
boundary hexagons in Fig. 1. 

We note that constraint operators commute not only 
with the full Hamiltonian but also with individual Q(abc)- 

[Pa,Q{abc)] = 0. The Hamiltonian (4) without constraint 
has an obvious ground state, |0), in which (t|Jj, = 1 for 
all rhombi. This ground state, however, violates the con- 
straint. This can be fixed noting that since operators 
Pa commute with the Hamiltonian, any state obtained 
from |0) by acting on it by Pa is also a ground state. 
We can now construct a true ground state satisfying the 
constraint by 

a ^ 

Here ^-^^ is a projector onto the subspace satisfying the 
constraint at site a and preserving the normalization. 

Obviously, the Hamiltonian (4) commutes with any 
product of Pa which is equal to to the product of cr^j, 
operators around a set of closed loops. These integrals 
of motion are fixed by the constraint. However, for a 
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topologically non-trivial system there appear a number 
of other integrals of motion. For a system with K open- 
ings a product of cr^^ operators along contour, 7 that 
begins at one opening and ends at another (or at the 
outer boundary, see Fig. 2) 

f, = J{<yl, (6) 

(7,) 

commutes with Hamiltonian and is not fixed by the con- 
straint. Physically these operators count the parity of 
'up' rhombi along such contour. The presence of these 
operators results in the degeneracy of the ground state. 
Note that multiplying such operator by an appropriate 
Pa gives similar operator defined on the shifted contour 
so all topologically equivalent contours give one new in- 
tegral of motion. Further, multiplying two operators 
defined along the contours beginning at the same (e.g. 
outer) boundary and ending in different openings, A, B 
is equivalent to the operator defined on the contour lead- 
ing from A to -B, so the independent operators can be 
defined (e.g.) by the set of contours that begin at one 
opening and ends at the outer boundary. The state \G) 
constructed above is not an eigenstate of these operators 
but this can be fixed defining 

iG/)=n^^^iG) (7) 

where Cq = ±1 is the eigenvalue of Tq operator defined 
on contour 7^. Equation (7) is the final expression for 
the ground state eigenf unctions. 

Construction of the excitations is similar to the con- 
struction of the ground state. First, one notices that 
since all operators Qabc commute with each other and 
with the constraints, any state of the system can be char- 
acterized by the eigenvalues [Qabc = ±1) of Qabc- The 
lowest excited state correspond to only one Qabc being 
— 1. Notice that a simple fiip of one rhombus (by opera- 
tor o'^ab) somewhere in the system changes the sign of two 
Qabc eigenvalues corresponding to two triangles to which 
it belongs. To change only one Qabc one needs to consider 
a continuous string of these flip operators starting from 
the boundary: \{abc)) = V(^abc)\^) with V(^abc) = D^' ^(cd) 
where the product is over all rhombi [cd) that belong to 
the path, 7', that begins at the boundary and ends at 
(a6c) (see Fig. 2 which shows one such path). This op- 
erator changes the sign of only one Qafcc, the one that 
corresponds to the 'last' triangle. This construction does 
not satisfy the constraint, so we have to apply the same 
'fix' as for the ground state construction above 

i-(«M)=n^^n^-('>Mio> (8) 

q ^ a ^ 

to get the final expression for the lowest energy excita- 
tions. The energy of each excitation is 2r. Note that a 
single flip excitation at a rhombus {ab) can be viewed as a 



combination of two elementary excitations located at the 
centers of the triangles to which rhombus {ab) belongs 
and has twice their energy. Generally, all excited states 
of the model (4) can be characterized as a number of 
elementary excitations (8), so they give exact quasipar- 
ticle basis. Note that creation of a quasiparticle at one 
boundary and moving it to another is equivalent to the 
Tq operator, so this process acts as t| in the space of the 
2^ degenerate ground states. As will be shown below, in 
the physical system of Josephson junctions these excita- 
tions carry charge 2e so that t| process is equivalent to 
the charge 2e transfer from one boundary to another. 

Consider now the matrix elements, 0^/3 = {Ga_\0\Gi3) 
of a local operator, O, between two ground states, e.g. of 
an operator that is composed of a small munber of (Tab- 
To evaluate this matrix element we flrst project a gen- 
eral operator onto the space that satisfies the constraint: 

6 VdV where "P = Jla The new (projected) 

operator is also local, it has the same matrix eleinents be- 
tween ground states but it commutes with all Pa- Since 
it is local it can be represented as a product of cr^ and 
Q operators which implies that it also commutes with all 
Tq. Thus, its matrix elements between different states 
are exactly zero. Further, using the fact that it com- 
mutes with Pa and Tq we write the difference between 
its diagonal elements evaluated between the states that 
differ by a parity over contour q as 

0+-0. = {0\\{^-^fqd\Q) (9) 

This equation can be viewed as a sum of products of az 

operators. Clearly to get a non-zero contribution each 
CT^ should enter even number of times. Each P contains 
a closed loop of six operators, so any product of these 
terms is also a collection of a closed loops of tr^. In con- 
trast to it, operator Tq contains a product of operators 
along the loop 7, so the product of them contains a string 
of operators along the contour that is topologically 
equivalent to 7. Thus, one gets a non-zero 0+ — 0_ only 
for the operators O which contain a string of a^' operators 
along the loop that is topologically equivalent to 7 which 
is impossible for a local operator. Thus, we conclude that 
for this model all non-diagonal matrix elements of a local 
operator are exactly zero while all diagonal are exactly 
equal. 



IV. EFFECT OF PHYSICAL PERTURBATIONS 

We now come back to the original physical system de- 
scribed by the Lagrangian (1) and derive the parameters 
of the model (4) and discuss the most important correc- 
tions to it and their effect. We begin with the derivation. 
In the limit of small charging energy the flip of three 
rhombi occurs by a virtual process in which the phase, 
(/>j at one (6-coordinated) island, i, changes by tt. In 
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the quasiclassical limit the phase differences on the indi- 
vidual junctions are <pind = ±7r/4; the leading quantum 
process changes the phase on one junction by 'Sn/2 and 
on others by —Tr/2 changing the phase across the rhom- 
bus (/) ^ <j) + TT. The phase differences, cj), satisfy the con- 
straint that the sum of them over the closed loops remain 
27r(r?, + <l>s/$o)- The simplest such process preserves the 
symmetry of the lattice, and changes simultaneously the 
phase differences on the three rhombi containing island 
i keeping all other phases constant. The action for such 
process is three times the action of elementary transitions 
of individual rhombi. So- 

r^Ey^E^Ui^pi-SSo), So = lAl^EjjEc (10) 

In the alternative process the phase differences between 
i and other islands change in turn, via high energy inter- 
mediate state in which one phase difference has changed 
while others remained close to their original values. The 
estimate for this action shows that it is larger than SS'o, 
so (10) gives the dominating contribution. There are in 
fact many processes that contribute to this transition: 
the phase of island i can change by ±7r and in addition 
in each rombus one can choose arbitrary the junction in 
which the phase changes by ±37r/2; the amplitude of all 
these processes should be added. This does not change 
the result qualitatively unless these amplitudes exactly 
cancel each other, which happens only if the charge of the 
island is exactly half-integer (because phase and charge 
are conjugate the amplitude difference of the processes 
that are different by I-k is exp(27ri(7) ). We assume that in 
a generic case this cancellation does not occur. External 
electrical fields (created by e.g. stray charges) might in- 
duce non-integer charges on each island which would lead 
to a randomness in the phase and amplitude of r. The 
phase of r can be eliminated by a proper gauge trans- 
formation I t)ob e*""''! X)ah and has no effect at all. 
The amplitude variations result in a position dependent 
quasiparticle energy. 

We now consider the corrections to the model (4). One 
important source of corrections is the difference of the ac- 
tual magnetic flux through each rhombus from the ideal 
value <i>o/2. If this difference is small it leads to the bias 
of 'up' versus 'down' states, their energy difference be- 
comes 2e = 27Ty/2^^Ej. Similarly, the difference of the 
actual flux through bavid's star and the difference in the 
Joscphson energies of individual contacts leads to the in- 
teraction between 'up' states: 

(ab) (ab),(cd) 

where Vab = e for uniform field deviating slightly from the 
ideal value and V(^ab).{cd) for rhombi belonging to the 
same hexagon. Consider now the effect of perturbations 
described by SHi, Eq. (11). These terms commute with 
the constraint but do not commute with the main term, 
H, so the ground state is no longer \G±). In other words, 
these terms create excitations (8) and give them kinetic 



energy. In the leading order of the perturbation the- 
ory the ground state becomes |G±) + ^ J2{ab) <^{ab)\^i±) 
Qualitatively, it corresponds to the appearance of virtual 
pairs of quasiparticles in the ground state. The density 
of these quasiparticles is |. As long as these quasipar- 
ticles do not form a topologically non-trivial string all 
previous conclusions remain valid. However, there is a 
non-zero amplitude to form such string - it is now ex- 
ponential in the system size. With exponential accu- 
racy this amplitude is (^)^ which leads to an energy 
splitting of the two ground state levels and the matrix 
elements of typical local operators of the same order 
E+-E_ ^0+~0- (^)^. 

The physical meaning of the W(ahc) excitations become 
more clear if one consider the effect of the addition of one 
operator to the end of the string defining the quasi- 
particle: it results in the charge transfer of 2e across this 
last rhombus. To prove this, note that the wave function 
of a superconductor corresponding to the state which is 
symmetric combination of | t) and | J.) is periodic with 
period tt and thus corresponds to charge which is multi- 
ple of 4e while the antisymmetric corresponds to charge 
(2n -I- l)2e. The action of one ctz induces the transition 
between these states and thus transfers the charge 2e. 
Thus, these excitations carry charge 2e. Note that con- 
tinuous degrees of freedom are charactc^rized by the long 
range order in cos(20) and thus correspond to the con- 
densation of pairs of Cooper pairs. In other words, this 
system superconducts with elementary charge 4e and has 
a gap, 2r, to the excitations carrying charge 2e. Similar 
pairing of Cooper pairs was shown to occur in a chain of 
rhombi in a recent paper^^; formation of a classical su- 
perconductive state with effective charge 6e in frustrated 
Kagome wire network was predicted in^^. 

The model (4) completely ignores the processes that 
violate the constraint at each hexagon. Such processes 
might violate the conservation of the topological invari- 
ants Tq and thus are important for a long time dynamics 
of the ground state manifold. In order to consider these 
processes we need to go back to the full description in- 
volving the continuous superconducting phases (pi. Since 
potential energy (2) is periodic in tt it is convenient to 
separate the degrees of freedom into continuous part (de- 
fined modulus tt) and discrete parts. Continuous parts 
have a long range order: < cos{2(f)o — 2(f)r) >'~ 1. The ele- 
mentary excitations of the continuous degrees of freedom 
are harmonic oscillations and vortices. The harmonic os- 
cillations interact with discrete degrees of freedom only 
through the local currents that they generate, further the 
potential (2) is very close to the quadratic, so we con- 
clude that they are practically decoupled from the rest 
of the system. In contrast to this vortices have an im- 
portant effect. By construction, the elementary vortex 
carries flux tt in this problem. Consider the structure of 
these vortices in a greater detail. The superconducting 
phase should change by or 27r when one moves around a 
closed loop. In a half vortex this is achieved if the gradual 
change by tt is compensated (or augmented) by a discrete 
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change by tt on a string of rhombi which costs no energy. 
Thus, from the viewpoint of discrete degrees of freedom 
the position of the vortex is the hexagon where constraint 
(3) is violated. The energy of the vortex is found from 
the usual arguments 

= ^ (In (i?) + c), c«1.2; (12) 

it is logarithmic in the vortex size, R. The process that 
changes the topological invariant, Tq is the one in which 
one half vortex completes a circle around an opening. 
The amplitude of such process is exponentially small: 
{i/Ey{D))^ where i is the amplitude to flip one rhom- 
bus and A is the length of the shortest path around the 
opening. In the quasiclassical limit the amplitude t can 
be estimated analogously to (10): ^EjEq exp(— 5*0). 
The half vortices would appear in a realistic system if 
the flux through each hexagon is systematically different 
from the ideal half integer value. The presence of free 
vortices destroys topological invariants, so a realistic sys- 
tem should either be not too large ( so that deviations of 
the total flux do not induce free vortices) or these vor- 
tices should be localized in prepared traps (e.g. David's 
stars with fluxes slightly larger or smaller than <3>s). If 
the absence of half vortices the model is equivalent to the 
Kitaev modeF placed on triangular lattice in the limit of 
the infinite energy of the excitation violating the con- 
straints. 

Quantitatively, the expression for the parameters of 
the model (4) become exact only if Ej ^ Ec- One ex- 
pects, however, that the qualitative conclusions remain 
the same and the formulas derived above provide good 
estimates of the scales even for Ej ^ Ec, provided that 
charging energy is not so large as to result in a phase tran- 
sition to a different phase. One expects this transition to 
occur at E* — rjEj with 77 ~ 1 which exact value can 
be reliably determined only from numerical simulations. 
Practically, since the perturbations induced by flux devi- 
ations from $0 are proportional to ^ and r becomes 
exponentially small at small Ec, the optimal choice of 
the parameters for the physical system is Ec ~ Ec- We 
show the schematics of the phase diagram in Fig. 3.'^^ 
The 'topological' phase is stable in a significant part of 
the phase diagram, further since the vortex excitations 
have logarithmic energy, we expect that this phase sur- 
vives at finite temperatures as well. In the thermody- 
namic limit, at T 7^ one gets a finite density of 2e 
carrying excitations {n^ ~ exp(— 2r/T)) but the vortices 
remain absent as long as temperature is below BKT-like 
depairing transition for half-vortices. 

V. QUANTUM MANIPULATIONS 

We now discuss the manipulation of the protected 
states formed in this system. First, we note that the 
topological invariant Tq has a simple physical meaning 
- it is the total phase difference (modulus 2tt) between 
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FIG. 3: Schematic of the phase diagram for half integer <E>s at 
low temperatures: S^d is the deviation of the magnetic flux 
through each rhombus from its ideal value. SC stands for 
usual superconducting phase, SCT for the phase with cos(2^) 
long range order of the continuous degrees of freedom and 
discrete topological order parameter discussed extensively in 
the bulk of the paper. The SCT phase and SC phase are 
separated by 2D quantum Ising phase transition. 

the inner and outer boundaries. This means that mea- 
suring this phase difference measures the state of the 
qubit. Also, introducing a weak coupling between these 
boundaries by a very weak Josephson circuit (character- 
ized by a small energy ej) would change the phase of 
these states in a controllable manner, e.g. in a unitary 
transformation = exp(iejir^). The transformation 
coupling two qubits can be obtained if one introduces a 
weak Josephson circuit that connects two different inner 
boundaries (corresponding to different qubits). Analo- 
gously, the virtual process involving half vortex motion 
around the opening gives the tunneling amplitude, et be- 
tween topological sectors, e.g. unitary transformation 
U'^ — exp(ztetT^). This tunneling can be controlled by 
magnetic field if the system is prepared with some num- 
ber of vortices that are pinned in the idle state in a special 
plaquette where the flux is integer. The slow (adiabatic) 
change of this flux to towards a normal (half-integer) 
value would release the vortex and result in the transi- 
tions between topological sectors with et ~ t/D2 

These operations are analogous to usual operations on 
a qubit and are prone to usual source of errors. This sys- 
tem, however, allows another type of operation that are 
naturally discrete. As we show above the transmission of 
the elementary quasiparticle across the system changes 
its state by r^. This implies that a discrete process of 
one pair transfer across the system is equivalent to the r| 
transformation. Similarly, a controlled process in which 
a vortex is moved around a hole results in a discrete Tq 
transformation. Moreover, this system allow one to make 
discrete transformations such as Vr^. Consider, for in- 
stance, a process in which by changing the total magnetic 
flux through the system one half vortex is placed in a cen- 
ter of the system shown in Fig. lb and then released. It 
can escape through the left or through the right bound- 
ary, in one case the state does not change, in another it 
changes by t^. The amplitudes add resulting in the oper- 
ation """^ . Analogously, using the electrostatic gate(s) 
to pump one charge 2e from one boundary to the island 



7 



in the center of the system and then releasing it results 
in a transformation. This type of processes allow 

a straightforward generalization for the array with many 
holes: there extra half vortex or charge should be placed 
at equal distances from the inner and outer boundaries. 

VI. PHYSICAL PROPERTIES OF SMALL 
ARRAYS 

Even without these applications for quantum computa- 
tion the physical properties of this array are remarkable: 
it exhibits a long range order in the square of the usual 
superconducting order parameter: (cos(2((^o — (f>r))) ~ 1 
without the usual order: (cos((/)o — 4>r)) = 0; the charge 
transferred through the system is quantized in the units 
of 4e. This can be tested in a interference experiment 
sketched in Fig. lb, as a function of external flux, $,„ 
the supcrcurrent through the loop should be periodic 
with half the usual period. This simpler array can be 
also used for a kind of 'spin-echo' experiment: applying 
two consecutive operations ^"^^ described above should 
give again a unique classical state while applying only 
one of them should result in a quantum superposition of 
two states with equal weight. 

The echo experiment can be used to measure the de- 
coherence time in this system. Generally one distin- 
guishes processes that flip the classical states and the 
ones that change their relative phases. In NMR lit- 
erature the former are referred to as transverse relax- 
ation and the latter as longitudinal one. The trans- 
verse relaxation occurs when a vortex is created and 
then moved around the opening by an external noise. 
Assuming a thermal noise, we estimate the rate of this 
process rj^ ~ f exp(— £'y(L)/T). Similarly, the transfer 
of a quasiparticle from the outer to the inner boundary 
changes the relative phase of the two states, leading to 
a longitudinal relaxation. This rate is proportinal to the 
density of quasiparticles, t^^ = TZexp{—2r/T). The co- 
efficient TZ depends on the details of the physical system. 
In an ideal system with some nonzero uniform value of 
e (defined above (11)) quasiparticles are delocalized and 
TZ ^ e/L2. Random deviations of fiuxes from half- 
integer value produce randomness in e, in which case one 
expects Anderson localization of quasiparticles due to off- 
diagonal disorder, with localization length of the order of 
lattice spacing, thus TZ ~ eexp(— cL) with c ^ 1, and e 
is the typical value of e. Stray charges induce random- 
ness in the values of r, i.e. add some diagonal disorder. 
When the random part of r, 5r becomes larger than e 
the localization becomes stronger: TZ ^ €{e/5r)^ where 
5r is the typical value of 6r. Upon a further increase of 
stray charge field there appear rare sites where Ti is much 
smaller than an average value. Such site acts as an addi- 
tional openings in the system. If the density of these sites 



is significant, the effective length that controls the deco- 
herence becomes the distance between these sites. For 
typical Ev{L) w Ej » 2K the transverse relaxation time 
reaches seconds for T ~ 0.1 iiT while realistic e/r ~ 0.1 
imply that due to a quasiparticle localization in a ran- 
dom case the longitudinal relaxation reaches the same 
scale for systems of size L ~ 10; note that temperature 
T has to be only somewhat lower than the excitation gap, 
2r, in order to make the longitudinal rate low. 

Most properties of the array are only weakly sensitive 
to the effect of stray charges: as discussed above, they 
result in a position dependent quasiparticle potential en- 
ergy which has very little effect because these quasipar- 
ticles had no kinetic energy and were localized anyway. 
A direct effect of stray charges on the topologically pro- 
tected subspace can be also physically described as a ef- 
fect of the electrostatic potential on the states with even 
and odd charges at the inner boundary; since the abso- 
lute value of the charge fluctuates strongly this effect is 
exponentially weak. 

Finally, we remark that the properties of the excita- 
tions and topological order parameter exhibited by this 
system are in many respects similar to the properties of 
the ring exchange and frustrated magnets models dis- 
cussed recently inii'i2>i9-29 



VII. CONCLUSION 

We have shown a Josephson junction array of a spe- 
cial type (shown in Fig. 1) has a degenerate ground state 
described by a topological order parameter. The man- 
ifold of these states is protected in the sense that local 
perturbations have exponentially weak effect on their rel- 
ative phases and transition amplitudes. The main build- 
ing block of the array is the rhombus which has two 
(almost) degenerate states, in the array discussed here 
these rhombi are assembled into hexagons but we expect 
that lattices in which these rhombi form other structures 
would have similar properties. However, the dynamics 
of these arrays is described by quartic (or higher) order 
spin exchange terms which have larger barrier in a quasi- 
classical regime implying that their parameter r is much 
smaller than in for the array considered here. This makes 
them more difficult to built in the interesting regime. 
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